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To  our  knowledge  the  results  of  the  pioneering  paper  on  the  boundary  layer 
theory  for  plates  by  K.  0.  Friedrichs  and  R.  F.  Dressier  have  not  been  extended 
to  anisotropic  plates.  Yet  the  basic  conclusion  of  the  above  paper,  namely  that 
the  equations  of  three  dimensional  isotropic  elasticity  can  be  split  into  two 
independent  systems  seems  to  have  been  accepted  by  some  authors  to  include  also 
the  case  of  anisotropic  elasticity.  While  this  happens  to  be  true  it  is  by  no 
means  obvious  as  the  following  presentation  illustrates;  hence,  we  believe  our 
providing  a  detailed  analysis  of  the  orthotropic  plate  will  be  of  interest  to  - 
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researchers  in  the  field  of  composite  materials  and  laminates. 
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Introduction 


To  our  knowledge  the  results  of  the  pioneering  paper  on 
the  boundary  layer  theory  for  plates  by  K.  0.  Friedrichs  and 
R.  F.  Dressier  El]  have  not  been  extended  to  anisotropic  plates. t 
Yet  the  basic  conclusion  of  the  above  paper,  namely  that  the 
equations  of  three  dimensional  isotropic  elasticity  can  be  split 
into  two  independent  systems  seems  to  have  been  accepted  by 
some  authors  to  include  also  the  case  of  anisotropic  elasticity 
T 7 ]  C  ].  While  this  happens  to  be  true  it  is  by  no  means  obvious 
as  the  following  presentation  illustrates;  hence,  we  believe  our 
providing  a  detailed  analysis  of  the  orthotropic  plate  will  be 
of  interest  to  researchers  in  the  field  of  composite  materials 
and  laminates. 

We  point  out  that  at  the  zeroth  order  of  approximation  of 
the  interior  solution  we  obtain  the  expected  thin,  orthotropic 
plate  equation  for  the  stress  function,  namely 


a22  a66  W'xxxx  +2  (2all  a2 2  '  2a12  "  a12  a66)  W,xxyy 


+  all  a66  w'yyyy  2  afi6  (an  a22  ~  ai2}  p(  )  (x'y) 


For  the  boundary  layer  approximation  in  the  limit  domain  we  find 
at  the  h^-step  that  the  stress  function  satisfies  an  equation 
having  a  three-dimensional  character  in  that  the  coefficients 
depend  on  the  Z-direction 


t  From  personal  correspondence  with  Professor  Friedrichs  we 
have  learned  that  he  also  is  unaware  of  any  generalizations 
of  this  type. 
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1  c  13  1  2 

TJ  (1_V23V32)  0  '  §  §  §  §  +  (-E^  +  V12V23J  6'§§§§ 


+  (1~vi2v21) 6'§§§§  "  0 


Using  typical  values  for  a  highly  anisotropic  material,  namely 

6  6  6 
Ej  =  20  x  10  psi,  E2  =  1.6  x  10  psi,  E^  =  1.3  x  10  psi , 

G12  =  °*6  x  1C)6  Psi'  G13  =  0,6  x  1q6  Psi»  G23  =  0,5  x  1q6  Psi* 
Vj^2  =  v^-j  =  v23  =  0.3  we  compute  the  above  coefficients  to  be 


a:  =  E ^  (1-v23v32^  =  0,7  x  10 


2b,  =  (-  4.  gi-  -  |-  V12V23)  -  1.6  X  10-6 


c:  =  E^  ^1-v12v21^  ~  0.046  x  10 


These  coefficients  show  that  the  equation  is  elliptic  as  the 

4  2 

roots  of  the  equation  an  +2bp  +c  =  0  are  seen  to  be  3  = 

+  i  (0.73)  and  v2  ^  =  +  i  (1.66).  In  this  case,  a  complex 
solution  may  be  sought  in  the  form  f 3 1  (pp.  29-30) 


F (x,y)  =  F-^x+^y)  +  F2(x+u2y)  +  F3  (x  )  +  F4(x+il2y) 


which  may  be  rewriteen  for  real  solutions  as 


F(x,y)  =  2  Re 


{f1(z1 


)  4  r2(z2) 


t  Provided  to  us  by  Professor  R.  B.  Pipes,  Director  of  the 
Composite  Materials  Center,  University  of  Delaware. 
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where  Z^:  =  x+p^y,  Z^’  =  x+i^y*  Consequently,  it  is  not 
technically  difficult  to  find  solutions  for  the  boundary  layer 
solution  using  function  theoretic  methods. 

Various  studies  on  interlaminate  stresses  have  been  made 
by  various  authors  [2] ,  [4-7] .  Most  of  these  essentially  avoid 

the  anisometric  nature  of  the  different  laminates.  It  would 
be  interesting  to  attempt  a  complete  analysis  of  the  laminate 
problem  bases  on  the  results  of  the  present  work. 


We  consider  in  this  work  the  case  of  an  orthotropic  elastic 
body*  namely  one  which  has  three  planes  of  elastic  symmetry.  Fur¬ 
thermore,  we  shall  assume  that  these  planes  are  mutually  perpen¬ 
dicular  and  that  they  are  orthoqonal  to  the  axes  of  a  Cartesian 
coordinate  system.  The  generalized  Hooke's  law  for  such  a  body 
may  be  written  as  f  3  ]  pg.  8. 


X  ' 

u'x  =  all°x 

+  a12°y 

+ 

al  3°  z 

(la) 

y:_ 

v'y  =  a12°x 

+  a22°y 

+ 

a23°z 

(lb) 

z 

w,z  =  a13°x 

+  a23ay 

+ 

a33°z 

(lc) 

yz: 

ii 

S 

+ 

N 

> 

II 

a,„T 

44  yz 

(Id) 

xz ' 

=  w'x  +  u'z  = 

9r  r  T 

55  xz 

(le) 

xy: 

=  U,y  +  v,x  = 

66  xy 

(If) 

The  elastic  constants  a^  may  be  expressed  in  terms  of  the 
engineering  constants,  (i=l,2,3)  (Young's  moduli), 

(i,j  =  l,2,3)  (Poisson's  ratio),  and  the  ^  (i,j=l,2,3)  (shear 
moduli)  as  follows: 
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‘12  “E.  '  a 
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Furthermore, 

Elv21  =  E2V12'  E2v 32  =  E3V23'  E3V13  =  Elv31 


(2) 


When  there  is  a  plane  of  isotropy  then  the  elastic  constants  simplify 
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The  problem  studied  by  Friedrichs  and  Dressier  [  ]  involved  an 

isotropic  body,  namely  the  case  of  complete  symmetry  where  E=E ' 
and  v=v ' .  In  what  follows  we  shall  show  that  their  analysis  for 
the  isotropic  plate  may  be  generalized  to  the  case  of  the  ortho¬ 
tropic  plate. 

We  consider  an  elastic  plate  of  uniform  thickness  2h,  lying 
parallel  to  one  plane  of  elastic  symmetry.  The  plate  shall  have 
an  edge  which  is  defined  by  a  continous  curve  B  possessing  a  con¬ 
tinuous  tangent.  The  weight  of  the  plate  and  other  body  forces 
are  to  be  neglected.  As  in  [  1  ]  we  consider  the  plate  to  be 
deformed  by  an  arbitrary  system  of  normal  stresses  distributed 
over  the  lateral  forces,  and  an  arbitrary  distribution  of  normal 
and  shear  stresses  which  vary  along  the  generators  and  perimeter. 
Adopting  the  notation  of  [  1  ]  we  consider  the  following  loads  to 
be  applied  to  the  plate 

a_(x,y,h)  =  -a(x,y)  top  face 

Z  (4) 

a  (x,y,-h)=  -b(x,y)  bottom  face 

Txz(x>V,±h)  =  Tyz(x,v,±h)  =  0. 

The  stresses  around  the  edge  are  given  by  o^(s,z),  ^ns(s'z)» 

t  _ (s, z) ,  where  t  (s,±h)=0  for  consistency. 
iiz  nz 

The  full  three-dimensional  system  of  equations  which  the 
plate  must  satisfy  are  given  by  , 
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+  T 

+  T  = 

X,X 

xy,y  XZ/Z 

+ 

a  +  r  = 

xy,x 

y,y  yz,z 

+ 

H 

+ 

Q 

II 

xz,x 

yz,y  z,z 

(5a) 

(5b) 

(5c) 


and  the  equations  (la-lf) .  The  equations  (la-lf)  involving  the 
displacements  may  be  replaced  by  six  compatability  equations 


-a.  .o  +2a, -o  +2a, -o  „_+(2a~, 
44  x,xx  13  x,yy  12  x,zz  23 

+2a33az/yy+(2a23+a44)oz/zz  =  0 

(2a13+a55)ax,xx+2all°x,zz+2a23ay,xx 
+2a33az,xx+(2a13+a55)oz,zz  =  0 

(2a12+a66)ox/xx+2all0x,yy+2a22oy,xx 
+2a13°z ,yy-a66az, zz  =  °' 

[2all°x+(2a12+a66) °y+ (2a13+a55}  az1  ' 


+a4  4) °y,yy+2a22°y,zz 

f 

-a55°y ,yy+2a12°y , zz 

i 

+(2a12+a66)ay,yy+2a23a 

yz+a44Tyz ,xx+a55Tyz ,yy 


(6a) 


(6b) 


z  ,xx 
(6c) 


+a66Tyz,zz 

[  (2a12+a66) °x+2a22ay+ (2a2  3+a44)  °z' 

+a^ r T  =  0  , 

66  xz,zz 

[  (2a13+a55)ax+(2a23+a44) ay+2a33az] 
+a66Txy,zz  “ 


( 6d) 

' xz+a44Txz ,xx+a55Txz ,yy 

(6e) 

'xy+a44Txy ,xx+a55Txy,yy 

( 6f ) 


The  thickness  of  the  plate  is  measured  in  the  z-direction. 
Following  [  1  ]  we  split  the  stresses  applied  around  the  edge  into 
even  and  odd  parts  with  respect  to  the  variable  z,  namely 


etc.;  similarly  the  face  conditions  are  rewritten 


-  -  (e) ,  -(0) 

o=o  +  a 
n  n  n 


in  terms  of  "even"  and 


o<0) (x,y,h) 
o|0) (x,y,-h) 
o^e) (x,y,lh) 


"odd"  functions  as 


P(x,v)  _  _a  b 

2  2  2  ' 


-  P(x,y)  _  a  b 

2  2  1' 


q(x,y)  _  _a  _  b 
2  "2  2  * 


(7a) 

(7b) 

(8) 


The  plate  problem  with  applied  boundary-forces  (s,z), 

x^(s,z),  x^(s,z)  with  face-forces  (7a, b)  will  be  called 
ns  n  z 

Problem  Ila  after  [  1  ] ,  whereas  the  plate  problem  with  applied 

.  /  p)  -  /  p  \  —  (  0  \ 

boundary-forces  o'  (s,z),  x  (s,z),  x'  (s,z)  and  plate 

n  ns  nz 

forces  (8)  will  be  called  Problem  lib.  The  Problem  Ila  is  a  pure 
bending  problem,  while  the  Problem  lib  is  a  generalized  plane 
stress  problem.  The  plate  problem  we  consider  may  be  uniquely 
decomposed  into  the  pair  of  problems  (IIa,IIb)  as  described  above. 

As  in  [1]  the  Problems  Ila,  and  lib  also  split;  indeed  we 
may  introduce  Problems  III,  IV,  V,  and  VI  such  that 

solution  III  +  solution  IV  =  solution  Ila, 

and 

solution  V  +  solution  VI  =  solution  lib. 

This  is  done  by  seeking  solutions  in  terms  of  even  and  odd  parts, 

6  0  0  0 

namely  u=u  +u  ,  a  =o  +a  ,  etc.  Upon  substituting  these  into 

XXX 

the  system  (5a-5c)  and  (la-lf)  the  resulting  equations  separate 
into  two  independent  systems,  namely  we  obtain 
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Problem  III 


0  0  0  0  0  0  ,,  . 
unknowns  f ax>  oy,  a z,  Txy,  u  ,  v  odd  m  z 


e  e  e 
Txy'  Tyz' 


even  m  z 


boundary  f (s,z),  Tns  (s»z)»  Tnz  (s,z) 
conditions 

-  o(0)  (x,y ,h)  =  -*jp,  a^0)  (x,v,-h)  =  ^p. 


Problem  IV 


unknowns  r  a®,  a®,  a®,  t®^,  u®,  ve  even  in  z 


0  0  0 
Txy'  Tyz' 


odd  in  z 


boundary 

conditions 


-e  n  -e  -0  _n 
°n'°'  Tns~° '  Tnz~° 


o®(x,y,lh)  =  0. 


Problem  V 


,  e  e  e  e  ,  e  e 

unknowns  {  °x'  °y '  °z'  xy'  U  '  v 

« 

0  0  o 

xy'  yz' 


even  m  z 


odd  in  z 


boundary  fa®(s,z),  t  (s,z),  t  (s,z) 


conditions 


o®(x,y,lh)  =  -y 


^Problem  IV  represents  a  rigid-body  motion. 


Problem  VI 
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.  000000  . 
unknowns  r  a  ,  o  ,  a „ ,  t  u  ,  v  odd  in  z 

I  x  *  y  z  xy 


e  e  .  o 
Txy '  T yz '  W 


even  in  z 


boundary  ,  -0  _  -0  -e 

conditions  f  °n~  '  Tns  '  Tnz  ' 


o°(x,y,±h)  =  0. 


Analysis  of  Problem  III 

Since  we  wish  to  obtain  asymptotic  results  as  h  >-0  we  introduce 
a  new  thickness  variable  c  =  z/h,  r,e[-l,l].  Likewise,  in  order 
that  stresses  and  displacements  will  not  vanish  or  become 


odd  in  t,  (9a) 

even  in  c ,  T^(s,il)=0.  (9b) 

nz 

On  the  top  and  bottom  of  the  plate  we  require  as  Friedrichs  and 
Dressier  [1]  that 

?Xz(x,y,ii,h)  =  Tyz (x,y,ll ,h)  =  0,  (10a) 

oz  =  -^p{3) (x,y)h3,  (10b) 

nz  =  ^-P(3)  (x,y)h3. 


unbounded  we  set  as  in  [  1  ] 


%  =  °n1} (s'r)h 


*ns  "  ^ns^5'^ 


TSz  =  Tnz)(s'c)h2 


(10c) 
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In  terms  of  the  x,y,t;  coordinates  our  system  (5a-5c)  ,  (la-lf) , 
(6a-6f)  become 


hg  +ht  +x  —  0> 

X,x  xy,y  lxz,c 

(11a) 

hT  +ho  +x  =  0, 

xy,x  y,y  yz,^ 

(lib) 

hx  +hx  +o  =  Of 

xz,x  yz,y  uz,^ 

(11c) 

U'x  =  allax+a120y+al 3°  z ' 

(12a) 

v'y  =  a120x+a220y+a2  30  z ' 

(12b) 

%  "  ha130x+ha230y+ha33°z ' 

(12c) 

ha44Tyz  =  v'C+hw'y' 

(12d) 

ha55Txy  =  hw'x+uV 

(12e) 

a66Txy  “  u,y+v'x' 

( 12f ) 

-a44h  ax,xx+2a13h2°x,yy+2a12ax,£;i;+(2a23+a44)h  °y,yy 

+2a22°y,??  +  2a33h2c,z,yy+(2a23+a44)oz,cc  =  °' 


(13a) 


(2a13+a55)h2°x,xx+2allax,cc+2a23h20y,xx-a55h  0y,yy 

(13b) 

+2a12°y,?c  +  2a33h2°z,xx'f  (2a13+a55)  °z,cc  =  °' 

<2a12+a66)h2°x,xx+2allh  °x,yy+2a22h2°y,xx 

+(2a12+a66)h2oy>yy+2a23h2oZf xx+2a13h2az>yy  (13c) 

"a66°z,^  = 


[  2a11hox+(2a12+a66)hoy+(2a13+a55)hozl ,yc+a44h  Tyzxx 


+a55^  Tyz ,yy+a66T vz , c c  ~ 


(13d) 


LI 


[(2a12+a66)hV2a22hV<2a23+a44)hoz,'xr. 

2 

+ac ,h  t  +a , , r  =  0, 

55  xz,yv  66  xz,rr 

(  (2a13+a55)h2ox+(2a23+a44)h2oy+2a33h2oz] 
2 

+a-ch  r  +ae/-T  -  0. 

55  xy  ,yy  66  xy  ,r.r. 


+  <l44^  Txz,xx 

( 1 3e) 


'xy+a44h  Txy,xx 
(13f ) 


Using  the  equations  (11,12,13)  we  seek  an  interior  expansion 
for  the  three  displacements  and  six  stresses.  We  assume  each  of 
these  unknowns  have  asymptotic  expansions  of  the  form 


u(x,y,c,h)  :=  £  u (x,y, c)hk,  etc. 

k=0 


(14) 


As  in  [  ]  we  obtain  for  the  hu  step  quite  directly  that 

,(0)  s  T(0)  s  o(0>  5  v(0)  s  u(0)s  T(0)  $  0  and  that  w(0)..M(x>y, 
xz  yz  z  xy 

is  independent  of  ^ .  The  equations  for  and  o become 

x  y 

all°x0)  +  a12°y0)  =  °' 

a12°i0)  +  a22°y0)  =  °' 

which  lead  to  c>  =  o  ^  =0  providing  that 
x  y 


_1_ 

'1E2 


aUa22'a12a2l  =  1717  a_V12v21)  *  °* 


The  analysis  for  the  h  -step  is  similar  to  that  in  [  )  so  we 

list  the  results 


r(1) 

5  T(1) 

S  a(1)  *  0 

xz 

yz 

z 

:=  W<1} 

(x,y) ,  v(1 

(1) 

T 

:=  T(1) 

(X,v)  . 

xy 

XV 

(1)  (  v  r  ii  ^  D  •  =11  ^  D 
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Furthermore,  providing  that  condition  (15)  holds  we  have  that 


(1) 


(1) 


:=  ^  (x»y)C.  :=  l  ( x , y ) c . 

x  v  y  v 


The  h  -step  leads  as  usual  to 


°i2)  E  °'  Txy>;*T^>(x'*,c'  Tyz>:“Ty2)(x'y)  (l~c2) 


=  T<2)  (x,y)  (l-c2)  ,  v(2)  :=V(2>  (x,y)C  , 


xz  xz 


u(2) :=U(2) (x,y) 5 , 


and  providing  conditions  (15)  is  met,  we  have  also 

m  (2)  (2)  (2) 

°x  =  I  (x,y)c,  a v  =  I  (x,y)c, 

x  x  y  y 


plus  a  new  result,  namely 

.(2)  .  1  (1) 


(1) 


w  =  2(al3  L(x'y)  +  a23  • 

x  y 


(0) 


(1) 


In  terms  of  the  function  W  it  is  possible  to  represent  £ 

.  <(2) 
and  l  as 


and 


(1) 


(1) 

l  (x,y) 

y 


- 1 - -  /  w(0)  +  w (o )  \ 

ana22_ai2  \  22  'xx  12,yy J 


- - 2  /r-a11w!3)  +  a,_W,(0)  ^ 

-a2  v  11  yy  12  xx / 

oo  an  N  / 


alla22  “12 


(15a) 


(15b) 


respectively.  Expressions  for  T*1*,  T^2*,  and  T^2^ 

xy  x  z  y  z 

be  found  in  terms  of  W^°'(x,y);  these  are 


may  also 
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(1) 


T =  — —  w|0)  , 

xy  a66  'xy 


T <*.y>  -  - 


2a66 (all“22  “12 J 


a22a66W'xxx 


+(2alla22“2a12“a12a66) 


W,(0)  *\ 
yyx  y 


(2) 


T^Mx,y)  - - - - —  ( 


2a  (a  a  -a2  )  l  alla66W'yyy 
2  66(alla22  a12  V 


(0) 


+<2alla22'2a12_a12a66)W 


(0)  \ 

'xxy  ) 


(15c) 


(15d) 


(15e) 


Using  these  in  (11c)  leads  to  the  following  differential  equation 
for  W(0) , 


a22a66W'xxxx+2(2alla22"2a12_a12a66)W'xxyy 


+a, i a- ,W 


lla66W'yyyy  =  “I  a^alla22~alb  P(3)^,y). 


( 1 5  f ) 


We  observe  that  equation  (15f)  is  the  equation  of  a  thin, 
orthotropic  plate  as  presented  in  the  book  of  Lekhnitskii  [  3  1 . 
It  is  interesting  to  note  that  this  reduces  to  the  result  of 
Friedrichs  and  Dressier  for  the  isotropic  plate. 


The  Boundary  Layer  Problem 

We  first  expand  the  edge  stresses  in  powers  of  the  arc  length 
variable s about  s=0,  and  then  replace  s  by  th 


°[1]  (s,Oh  =  o  (0 ,  c) h  +  (0  r r. )  th2  +  .. 


n ,  s 


'ns,(s’')h  -  ^a)(0-tlh  +  ’naU(0',;>th2  *  " 


(16a) 

(16b) 


riV  <s,Oh2  =  7<2)  (0,C)h2  +  tJ[2)  (0,c)th3  +  ... 


nz 


nz 


P*  _ 


(16c) 
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and  the  face  stresses  as 

°z  =  ^X'Y'-l'h)  =  +P(o!o)h3  +  *•*  '  (17a) 

Txz  "  Tyz  "  0  for  «  "  il-  (17b) 


If  the  boundary  curve  B  is  parametrized  by  x=x(s),  y=y(s)  then 

under  the  transformation  £=x/h,  n=y/h  we  obtain  a  new  boundary 

curve  B*  with  arc  length  t=s/h,  where  the  parametrization  is 

now  given  by  c=£*(t,h),  n=n*(t,h).  The  unit  tangent  and  normal 

vectors  are  given  by  [X  (t,h) ,Y  ]  ,  [X_.Y  ]  respectively.  In 

s  s  n  n 

this  notation  the  edge  conditions  may  be  rewritten  (see 
Lekhnitskii,  pg.  2),  as 

01!1)  (0,Oh  +  o*1*  (0,t)th2  +  ... 
n  n  f  s 

-  x^>*'"*'‘-hl+2xnyntxy*VnV  (18a 

<°«>h  +  ’ns!s(0'',h2t  +  ••• 


=  X  Y  o  +  [X  Y  +X  Y  ]  t  +Y  Y  a  , 
nszlns  s  n’  xy  nsy' 


;{„lUa,Oh2  *  t<2)rS(o,c)th2  +  ... 


=  Xn*xz  +  Yr/yz 


(18b) 


(18c) 


Here  we  have  used  the  notation  of  [  ]  to  indicate  that 


o(f.  rh) 


o (he , hn » c  »h) ,  etc. 


In  what  follows  we  match  interior  and  exterior  expansions  to  obtain 
the  asymptotic  behavior  of  the  various  terms  in  the  plus  super¬ 
script  function  expansions  as  {>-»;  for  details  of  the 


procedure  the  reader  should  consult  [  1  1 


The  limit  domain  in  the  £,n,C  space  is  given  by  D*:={(£,n,c) 


E.<0,  |n|<“,  -1<_c<_1).  We  illustrate  this  by  showing  its  cross- 
section  in  the  £-c  plane  as 


The  segment  AB  corresponds  to  the  boundary  edge,  whereas  BC  and 
AD  correspond  to  the  top  and  bottom  faces  respectively. 


The  h°  step 

In  the  limit  domain  in  the  £,n,C -space  we  obtain  the 
following  boundary  conditions  at  the  two  faces  and  the  edge 


+  (0)  ... 

0  _  — 


+  (0) 
Tyz 


+  (0) 
xz 


=  0 


+  (0)  _ 
0  __  ““ 


+  (0) 
Txy 


_  +(0)  _ 


xz 


=  0 


on  5  =  1 1 
on  £  =  0 


(19a) 

(19b) 


At  £=-«  o ^  (£ , n f C ) -o  ^  (0 * 0 , £ ) .  Furthermore,  since  for  the 

interior  problem  h0  we  obtain  the  boundary 

x  x  y  x  z 

conditions 


+  (0)  _ 
ffx 


-  r 


+  (0)  _  +(0)  _ 


xy 


-  r 


xz 


=  0  at  £  =  -o 


(20) 


As  in  [  ]  one  has  the  conditions 
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1 


9n 


0, 


32S(1) 


=  0, 


3no (n-1) 


3n 


n 


5  o. 


(21) 


which  leads  to  a  splitting  of  the  stress  system  into  a  plane 

strain  problem  involving  six  equations  in  the  four  unknowns 

'  K°]  '  an^  *-nto  a  torsion  problem  having  three 

equations  in  the  two  unknowns  ,  t  ^ .  Obtaining  the  asymp- 

x>  yz 

totic  conditions  for  0^®^  ,  from  the  interior  solu¬ 

tions  and  integrating  the  zero*"*1  under  system  yields 


<*  *(0)+a  t{0) 

a12°x  a23°z 


)  • 


Since  all  applied  stresses  are  zero  on  the  boundary  of  the 

limit  domain,  we  conclude  as  in  [1]  that  a*0^,  '  are 

identically  zero.  This  implies  that  ciy0^  =  0  also.  The  torsion 

problem  has  vanishing  boundary  conditions  and  therefore  it  follows 

that  t(0),t(0)-r0. 
xy  yz 


The  h*  Step 

We  consider  next  the  plane  strain  problem  for  ,t^  , 

X  X  z 

t'U  ,0^^ •  Since  all  n  derivatives  vanish  we  are  led  to  the 
following  six  equations  (Here  we  have  suppressed  the  +(1)  super 
scripts  in  our  notation.): 


°x,s  +  Txz , C 
Txz,C  +  °z,q 


0, 

0, 


(22a) 

(22b) 


"a44ax,5  5  +  2a12°x,j;J;  +  2a220y,j;i;+(2a23+a44)oz,t;c 


=0,  (22c) 


1  7 


(2ai3+a55)ox,£{;+2al10x  ,CQ  +  2a23aytF,r, 
+(2a13+a55)oz ,cc  =  °' 

(2a12+a66)ox,^  +  2a22uy,a  +  2a230z,eC 

I(2a12+a66)ox+2a22ay+(2a23+a44)oz]' 

+a66Txz,cc  =  0  ' 


+2al 2°v , 


r  r 


+  2a 


"a66l,z,cc"°' 

SC+a44Txz,£<- 


33°z,CS 

( 22d) 


( 22e) 
( 22f ) 


Equations  (22a, b)  imply  the  existence  of  a  stress  function 

0:=©(C,O  such  that  ax=:0'tz'  °z=:0's£'  Txz=:-0'EC*  From 

this  eouations  (22c-22f)  may  be  rewritten  as 


■a440'cas  +  2a120'c^c  +  2a22°v,cc  +  (2a23+a44)  0,f:m  °' 

(2ai3+a55,0',a^2a11G'c^c+2a23oy^e+2ai2°y,^ 

+2a330'CCCC+ (2a13+a55) 0 ' ~  °' 

(2a12+a66)0,?c5E  +  2a22oyf  c?+2a230,s??c-a660,5ec!.  -0, 
[  (2a12+a66)0,??  +  2a22ay+(2a23+a44)0,^J 

'0,Ca?a44“a660,CCCi;  =  °' 


(23c) 


(23d) 

(23e) 


( 2  3  f ) 


Integration  of  (23f)  leads  to 

a22°y  =  -a120'cra230'CC  +  *l(C)+*2U) 

=  -a120x~a23°z+hU)+*2U)  ' 

Using  the  interior  expansions  and  matching,  one  obtains  the 
asymptotic  behavior  as  £•*—«» 


from  which  follows 
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a22° y  =  -a12°x  "  a23°z  +  cl?  +  c2?  +  c3;  (2 

however,  as  is  odd  in  the  c -variable  we  are  forced  to  have 

=  c3  =  0.  After  some  further  manipulation  with  the  equations 
(23c-e)  it  is  found  that  0  satisfies 


(a 


2 

22a33-a23,0'£££S+(a22 
+  (alla22~a12^  0  '  f;  £  £  £ 


[2a13+a55] -2a12a23) 0 
=  0. 


(25) 


For  the  case  of  an  isotropic  material  equation  (25)  reduces  to 
the  biharmonic  eauation  as  was  shown  in  the  paper  of  Friedricks 
and  Dressier.  For  a  material  with  one  plane  of  isotropy  the 
equation  becomes 


<i-§T*l2>e.Ee£C  +  2U-vV)0,tuc  ♦  ^d-v2)e,„„-o.  (26) 

whereas  for  an  orthotropic  material  (25)  mav  be  expressed  in 
terms  of  the  elastic  coefficients  as 

1  2v13  1  2 

EJ(1_v23v32)0'££££+(  ITj-  +  “  E^  V12V23)  0,££CC 

,  (27) 

+  EJ(1-V12V21,0'CCCC  =  °* 

Equations  (26,27)  show  the  three-dimensional  nature  of  the  equation 
which  0  satisfies  in  that  it  contains  elastic  constants 
associated  with  the  z-direction. 

Recalling  that  cr  *:0,rr,  a  =:0,fr,  -xv  =:0,r,  we  obtain 

A  s s  Z  ss  XZ  ss 

the  following  boundary  conditions  for  0,  namely 


10 


°'CC  = 


°’U  -  0 


on  AB 


(28a) 


and 


0,  =  0,?£.  =  0  on  BC  and  DA 


(28b) 


On  CD  we  obtain 


0, 


(1) 

l  (0,0)  C,  and  0,^  =  0,^  =  0. 


(1) 


From  (24),  (29)  and  the  fact  that  a  =  7  (0,0)^  on  CD  the 

y  L 

x 

expression  for  a  becomes 


.(1) 


(1) 


(29) 


l22ay  a120,^c_a230'Cc+(a22  ^  (0,0)+a12  \^0,0))z,. 


This  completes  the  h  -step  for  a  ,a  ,t  ,  and  o  ,  modulo  the 

x  z  x  z  y 

fact  that  our  plane-strain  problem  is  defined  only  in  terms  of 
the  "interior"  stress  problem  evaluated  at  the  edge  of  the  plate. 
For  further  details  on  how  to  determine  the  proper  boundary  con¬ 
ditions  for  the  interior  problem  in  order  to  determine 

U)  ,(1) 

I  (0,0),  l  (0,0)  the  reader  is  referred  to  [1],  pp.  21-22. 
x  x 

We  turn  now  to  the  torsion  problem  for  '  ^yz^ '  as 

above  we  shall  omit  the  +(1)  superscripts. 

Equations  (lib) ,  (13d) ,  (13f)  reduce  to 


Txy , C  +  Tyz,c  °' 


a44Tyz  ,  c  £  +  a66Tyz ,  x,  c  =  °' 


(30a) 

(30b) 


a44Txy'U  +  a66Txy,^  =  °' 


(30c) 
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Introducing 


T 


xy 


-T 


yz 


$ 


'£ 


allows  (30b, c)  to  be  rewritten  in  the  form 


(31) 


a44^  '  £  £  £  +  a66  ^'££5  “ 
a444’ '  ^  CC  +  a66*'??i ;  =  °' 
which  imply  that 

a44*'££  +  a66 =  °' 


(32a) 

(32b) 


(33a) 


where  6  is  an  arbitrary  constant.  This  reduces  to  the 
Friedricks-Dressler  form,  A4>=<5,  when  the  material  is  isotropic. 
When  there  is  a  plane  of  isotropy  this  becomes 


E(l+v')*,cc  +  E'(l+v)*  =  6  ,  (33b) 

and  for  an  orthogropic  material  equation  (33a)  in  terms  of  its 
elastic  coefficients  is 


5  . 


The  boundary  conditions  for  (33a)  are 

♦  s-  Txy)  ( 0 » n » C )  =  (0,C)  on  AB  (34) 

and 

4>,c  :=  -Tyg*  U » n  r-1 ,h)  =0  on  BC  and  AD.  (35) 

Since  -*■  T^^  (0,0) q  as  £-*■-»,  we  have 

♦  »_  -  T***  (0,0)  6  on  CD.  (36) 

c,  xy 

This  boundary  value  problem  may  be  solved  by  Fourier  series  as 
suqgostod  in  f  1  1.  The  arbitrary  coefficient  is  determined  to 
be  T*1,)  (0,0). 


21 


The  h2-Step 

Here  the  boundary  conditions  are,  as  usual,  the  same  as  in 

[  ];  the  reader  is  referred  to  equations  (51-55)  of  this  work 

for  details.  We  treat  next  the  case  of  plane-strain  for  the 

so-called  "excess"  stresses,  [1],  pg.  25,  o  ^ and 
+  ( 2 ) 

a  .  As  has  been  our  custom  we  introduce  a  stress  function 
—  z,  n 

r  (£  ,  x, )  such  that 


r'?^:  — x, n ' 


r'cr"-z,n' 


r , 


U 


:  =— t 

-xz  ,y 


(37) 


Then  r  is  seen  to  satisfy 


(a 


2 

22a33~a23) r,^CCC+(a22 
+  (alla22'a12) F  'CCCC 


[2a13+a55] -?a12a23) r 
=  0. 


cm 


(38) 


Furthermore,  a  is  found  to  be 

-y,  n 


l22— y ,  n  'a12-x,n"a23-z,n+4,l(C)+<t,2(c) 

=  ~a12r  'qc~a23r  '££+<(,l +<t>2  5 


(39) 


From  the  fact  that  «  is  odd  in  c  and  all  the  stresses 

— y  »n 

vanish  at  infinity  we  conclude  ^1-^2=° - 

+  ( 2. ) 

We  turn  next  to  the  "excess"  torsion  problem  for  t '  '  and 

xy  f  n 

+  ( 2 ) 

Lyx  t)  an<*  f°How^n9  our  previous  notational  abreviation  drop  the 
+(2)  superscript.  We  have  from  (30a-c) 


Txv,nC  +  Tyz,n  r, 
a44Tyz,nC?  +  a66Tyz,nc?  _ 
a44Txv,nCC  +  a66Txv ,ncc  =  °* 


(40a) 

(40b) 

(40c) 


i  if?*  •»*■<■' 
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Setting  <t>  ,.  :=  *  $  ,r  :=  -t  ,  we  obtain  after  integration 

r  C  xv,n  'C  yz  ,n 

and  considering  the  behavior  as  f>-«,  that 


+  a 


66  'CC 


=  0. 


(41) 


Again  the  "excess"  stress  can  be  determined  by  solving  (41)  with 
<p  prescribed  over  the  boundary  of  our  fundamental  strip.  We 
must  now  consider  the  system  of  "excess"  stresses  themselves, 
namely  functions  of  the  form 


J.(2)(^ru?)  :=  S,(2)  (£,5>n  +  A.  (e,n)  .  (42) 

,  n  1 

The  eauations  for  the  A^(£,n)  are  obtained  by  direct  substitution 
into  the  nine  elasticity  equations  in  the  e;  ,  n ,  c  defining  the 
(+) 2  excess  stresses.  As  before,  the  system  splits  into  two 
independent  sets  one  which  we  refer  to,  after  Friedrichs-Dressler , 
as  "quasi-plane-strain"  and  the  other  as  "quasi-torsion." 


Quasi-plane-strain 

From  the  equations  (11a) ,  (11c) ,  (13a) ,  (13b) ,  (13c) ,  and 
(13e)  we  obtain 

+  ( 2) 

Ax.C  +  flxz,c  “  -Txy,„  -*’s>  (43a) 

*xz,e  *  Vs  -  4yz!n  -=  <43b) 

“344Ax,  f;£  +  2a12Ax,  Cj;  +  2a22Ay,  et+(2a23+a44)  Az,  CC  =  °»  (43c) 


(2a13+a55)  AX,£5  +  2allAx,  ?i;  +  2a2  3Ay,  5^  +  2a12Ay,  Ci; 
+2a33Az,C5+(2a13+a55)Az,c^  =  °» 


(43d) 


(2a, _+a,c) A„  „+2a,,0A,,  __+2a~,A„  „  -a^A,.  =  0, 

12  66'  x,^f;  22  y,£,Z  23  z,CC  66  z,i;t  ' 


(43e) 
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(  (2a12+a66)  Ax+2a22Ay+  (2a2  3+a44)  Az^  '  F,F+a4  4Axz  ,  FF 
+a66Axz , f  f  =  °* 


For  an  orthotropic  material  these  simplify  somewhat,  namelv 


(4  3f ) 


equations  (43c-f)  become 


( 1  +  v  ' )  A  +  ( 1  +  v )  E—  A  -~(A  +  A  ) -A  =0, 

x  ,  F  F  E  x  ,  f  r,  E  x  ,  r,  f  v  ,  f  F  z  ,  r,  f  ' 


( 4  4c) 


<ltv,,Vu+atv,jr 


Az  ,  c  F  ° ' 


F.uT'VVr  V.«  *  •< 


( 4  4d) 


( 4  4e) 


<VVr-  *  “•  {44f) 

which  in  turn  reduce  to  that  of  [  1  ]  for  the  isotropic  case.  The 
boundary  conditions  on  the  ciuas i-plane-strain  are  the  same  as 
those  in  Friedrichs-Dressler ,  pp.  28-29. 


Quasi-torsion 

From  eauations  (lib),  (13d),  (13f)  we  obtain 


Axy,f.  +  Av  z ,  r.  -0v,n  a22  f'ff  +  a22  r  '  F  F  ' 


a44Avz,FF+a66Ayz,r.F  (2all'a22  <2a12+a66)  r,FFF 


+  (2a13+a55“a22(2a12+a66) )  r  *  F F F 1  ' 


(45a) 


(45b) 
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a44Axy,U+a66Axy,ee  =  "  f  (2a13+a55-a22  (2a23+a44)  }  r  'ccC 


a 


(45c) 
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+  (2a33-a22(2a23+a4  4)  )  r  '  c  C  f;  * 


Eauations  (45b,c)  reduce  to  the  Friedrichs-Dressler  case  for 
isotropic  materials,  namely  we  have 


and 


AA 


yz 


-AT , ? , 


(45b') 


AA 


xy 


(45c’) 


d 
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